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Abstract 

Using the recent formula presented in [1, 2] for the hnk representation of tree- level Af = 8 su- 
pergravity amplitudes, we derived a CSW-like expansion for the Next-to-MHV 6- and 7-graviton 
amplitudes by using the global residue theorem; a technique introduced originally for Yang-Mills am- 
psj plitudes [3]. We analytically checked the equivalence of one of the CSW terms and its corresponding 
Risager's diagram [4]. For the remaining 6-graviton and all 7-graviton terms, we numerically checked 
the agreement with Risager's expansion. We showed that the conditions for the absence of contribu- 
fSJ tions at infinity of the global residue theorem are satisfied for any number of particles. This means 
^ that our technique and Risager's should disagree starting at twelve particles where Risager's method 
CN is known to fail. 
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1 Introduction and Main Results 



In their seminal work [5], Cachazo, Svrcek and Witten (CSW) proposed an analytic expansion for 
tree-level Yang-Mills amplitudes. The rules of the expansion are similar to Feynman rules except for 
the vertices which are off-shell continuations of Maximally- Helicity- Violating (MHV) amplitudes. The 
off-shellness is provided by introducing an auxiliary spinor into the calculations. Risager [4], inspired 
by the Britto-Cachazo- Feng- Witten (BCFW) recursive method [6, 7], later presented a proof of the 
CSW rules to compute gluon amplitudes. In his approach, an auxiliary spinor similar to CSW's is 
used to apply a certain complex deformation on external momenta. Like in the BCFW method, the 
residue theorem is used to recover the non-deformed Yang-Mills amplitude. 

The simplicity of the CSW expressions in gauge theory motivated similar computations for gravity. 
However, when applied to gravity, Risager's method only produces graviton amplitudes up to eleven 
particles in the Next-to-MHV sector. The failure of Risager's method to write an MHV-expansion 
for NMHV gravity was first discovered by Bianchi et al. [8] using numerical techniques and later 
analytically confirmed by Benincasa et al. [9]. Also, recently the difference of the pure 12-graviton 
NMHV amplitude with its Risager's expansion was analytically calculated by Conde and one of the 
authors of this paper [10]. 

Here we briefly explain Risager's method. Throughout this paper, we only concentrate on pure 
(gluon or) graviton amplitudes in the NMHV sector which contain three negative helicity and n — 3 
positive helicity particles. Using the spinor-helicity notation for the on-shell momentum of particle 
i, = A* A*, Risager's expansion is obtained from a complex deformation on the anti-holomorphic 
spinors. A*, of the external negative helicity particles. Since we restricted ourselves to the NMHV 
amplitudes, Risager's deformation with a complex variable z is given by 

AW(z) = Xi") +z{bc)fi 

AW(z) = AW + z(ca)7/ (1.1) 
AW(z) = A(^) + z (06)77 

where the negative helicities are labeled by a, b and c, and fj is an arbitrary spinor. 

Through this deformation, the amplitude becomes a rational function of z, Ainiz)- One can apply 
the residue theorem to construct the amplitude Ain from the poles of M.n{z). Risager's expansion, 
J^n^^, is obtained by summing the residues of M-n{z): 

Mf^=Y,ML (a-,L+,(-p)') ^^^^^^^^^, A^^ (p+,6-,c-,i?+) , (1.2) 

where P = pa + Pl+ is Risager's deformed momentum flowing in the internal propagator, and 
{R^) denotes the subset of external positive- helicity particles in the left (right) sub- amplitudes in 
(1.2). The hatted momenta are evaluated at the position of the poles which makes the corresponding 
propagator on-shell. Also, note that each term in the expansion is the product of two MHV amplitudes 
and a propagator. 

As a consequence of the residue theorem, M.ri{z) must vanish at infinity for the expansion to 
be valid, but the behavior for gravity is in fact Mn{z) ~ z"^"*^^ when z — )■ 00. This MHV-vertex 
expansion for gravity, hence, is not valid when applied to amplitudes of more than eleven gravitons. 

As shown in [11], Yang-Mills MHV amplitudes are localized on lines in Z supertwistor space. The 
CSW expansion in Yang-Mills tells us that NMHV amplitudes are localized on pairs of lines in twistor 
space. Gravity MHV amplitudes, on the other hand, are not localized on lines because they have 
derivative of a (5-function support on lines. However, in the NMHV sector one can write amplitudes 
as integrals over the Grassmannian G(3,n) which is the space of 3-planes in C". A convenient way 
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to represent a point in G(3, n) is as a 3 x n matrix modulo a GL(3) action. This matrix can be 
thought of as n 3-vectors which are the columns of the matrix. The columns are denoted by Ca G C^, 
or ec[uivalently [c^^ 1 ~ as they are non-zero. We call the space of 3-vectors the C-space and will 
discuss localizations in this space throughout this paper. 

In Yang-Mills, NMHV amplitudes are also localized on pairs of lines in the C-space. Hence, lo- 
calization in twistor space and CP^ are equivalent for Yang-Mills. We take CSW in gravity to mean 
that NMHV amplitudes are localized only on pairs of lines in the C-space. For this reason, we call 
the formulation in this paper a CSW-like expansion for gravity. 

Until recently, there has been an obvious lack of tools to explore a CSW-like expansion in gravity. 
One problem has been the absence of compact expressions for MHV amplitudes, the building blocks. 
This issue was overcome when Hodges presented an elegant formula for the MHV gravity amplitudes, 
analogous to the Parke- Taylor formula for Yang-Mills [12]. Another obstacle has been the lack of a 
Grassmannian formulation. 

In [13], Cachazo and Skinner proposed a manifestly permutation invariant formula to compute all 
tree- level = 8 supergravity amplitudes in terms of higher degree rational maps to twistor space. The 
parity invariance and soft limits of the formula, two strong evidences for the proposal to be correct, 
were checked in [14]. This conjecture was later on proved in [2], showing that the formula admits the 
BCFW recursion relations and also produces 3-particle MHV and MHV amplitudes. Moreover, the 
formula behaves at large momenta (large BCFW's complex variable, z) as tree-level gravity amplitudes 
do, l/z^ [15]. In [1] and independently in the same work of Cachazo-Mason-Skinner [2] the formula 
was presented in the link representation. This formulation was indeed another missing ingredient to 
further explore CSW in gravity. Cachazo-Skinner formula and the link representation are reviewed in 
section 2. 

In this note, we are looking for a formula analogous to the CSW one but for gravity. The analogy 
is not directly in twistor space but is in C-space. The main procedure is explained in section 3. 
We considered the link representation of the Cachazo-Skinner formula and applied a global residue 
theorem (CRT) to it [16]. The residue theorem writes the amplitude as a sum over terms which 
coincide with different localizations of particles on a pair of lines in the C-space. These are in fact 
the CSW terms for NMHV gravity amplitudes. 

We computed the CSW terms for the NMHV 6- and 7-graviton scattering amplitudes in sections 
4 and 5, respectively. The results obtained are compared to the known corresponding Risager's 
diagrams. We show that the two sides are in complete agreement. Both theories, i.e., Yang-Mills and 
gravity, have an expansion in terms of two lines in C-space. Throughout this paper, we graphically 
represent each term in (1.2) for gravity as two intersecting lines as shown in figure 1. 

h- 

Figure 1: The CSW-like localization of NMHV amplitudes on lines in C-space. 

For 6- and 7-graviton amplitudes, there are 21 and 45 Risager's diagrams respectively, which are 
the same as the localizations of the Cachazo-Skinner formula. This formula however, does not hold 
us back from exceeding 12-graviton amplitudes. Indeed, a power counting in section 6 shows that the 
global residue theorem is valid for any number of particles. In section 7 we discuss a possible reason 
for the discrepancy between Risager's and our technique. 
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2 Twistor String and Link Representation 

Recently, a new formula for the tree-level S-matrix of = 8 supergravity for all R-charge sectors 
(labeled by d = + 1) was proposed by Cachazo and Skinner [13]: 

-Mn = Vor(GL(2))" ^^^^ ^^^^ n '^'^^ ^'(^^ - ^(^i))e^P ([/^(^O A.] +xV) , (2.1) 

whose ingredients are explained below: 

• ^ is a holomorphic map from an n-punctured Riemann sphere S with homogeneous coordinates 
{a-, CJ-) to = 8 supertwistor space CP^'^ given by degree d polynomials 

d 

Z'{a) = Y,2i{a'-r{a'-f-'', (2.2) 

a=0 



$ and $ are n x n matrices defined as 



and 



n,^j(ifc)i>„(»j'") 



j^t a=0 

where pa are reference points in S. As written above, ^ has rank d while $ has rank d = n—d—2. 
In order to obtain the reduced matrices ^^'^'^ and ^^'^'^ with non- vanishing determinant, one must 
first delete n — d rows and columns of and d + 2 rows and columns of <I>. After this, the two 
det' in (2.1) are defined as: 

Icbredl _ |$red| 

det'(«>) = i ^ r, det'($) = ^ ^ r, (2.5) 

In . . . rrfllci ...Cd\ In . . . rd+2| |ci . . . Crf+2| 

where |ri . . .r^l and |ci . . .Crf| are the Vandermonde determinants of the remaining rows and 
columns for $ and |ri . . . rdj^2\ and |ci . . . Cd+2\ of the deleted rows and columns for $: 



|ci...q|= n {ij), (2.6) 

i j' £ {remaining} 

4: |ci...Q+2|= n (^J')' (2.7) 
and the same for |ri . . . r^l, |ri . . . r(i+2|- 



i,jS{deleted} 



4 



It is known that (2.1) is completely localized by the 5-functions, being the analog of the twistor string 
formulation for = 4 super Yang-Mills tree amplitudes proposed by Witten [11] and later studied 
by Roiban, Spradlin and Volovich [17]. However, the constraints involve solving polynomials of high 
degree and summing the contribution of all solutions. 

To simplify this problem, He [1], Cachazo, Mason and Skinner [2] wrote the amplitude in the 
so-called link representation, which is a "gauge fixed" Grassmannian G{k, n) formulation in which 
k = d+ 1 of the external particles are taken to be supertwistor space Z eigenstates while the remaining 
ones are taken to be dual supertwistor space W eigenstates. For amplitudes having only gravitons as 
incoming particles, we choose the gauge such that the k columns associated with the negative helicity 
gravitons form an identity matrix. We label the gauge fixed columns by r and the remaining ones by 
a. The unfixed variables Cra are called link variables. 

In summary, the wavefunctions of the particles labeled by r are localized in Z space, while the 
wavefunctions of the particles labeled by a are localized in W space. In the link representation, each 
3x3 minor becomes linear in the link variables, and the price to pay is that the amplitude is now a 
contour integral in {d — l){d — 1) variables. 

As presented in [2], the final product for the N'^^^MHV tree amplitude in the link representation 
reads 



Mn,d iZr,Wa) 



n 



dc 



ra 



D\ 



n—ln 
12 



naClaC2anr 

C-rn—l^r 



W ^12 Ci(jC2aCrn— iCf 

{ab) 



Tjab 



Trn—ln 
Iirs 



^yl2r 
' an— In 



(2. 



exp (iCraZ'^W) , 



where 



jjab 

rs 



Csa 



Csb 



ab 



D 



ab 



ValLin = (12r)(ran - l)(n - lnl)(2on) - (2 r a)(a n - 1 n)(n 1 2)(rn - 1 1) 

We call V^^1i„ the Veronese polynomial and the reason for the name will be explained in the next 
section. (/)^'^\ (p^^^ are the determinants present in (2.1) with the following particular choice of removed 
rows and columns: using cyclicity, the labels are chosen such that particle 1 has negative helicity and 
particle n has positive helicity. Then, 

• from (2.3), all rows and columns a plus r = 1 are removed, 

• from (2.4), all rows and columns r plus a = n are removed. 

In other words, ch^ ' is the determinant of a d x d matrix with elements 



C-raCrb^saC-sb 

{2r a){an — 1 n)(n 1 2)(rn 



(2.9) 
(2.10) 



Trn—ln 
ilrs 



for r 7^ s and r, s 7^ 1, 



/ -J Trn—ln ' 
ilr 



(2.11) 
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where as before, r and s run over the particles labeled by Z leaving out 1, and (j)^ ' is the determinant 
of a d X d matrix with elements 

^"■^ = !^ for a / 6 and a,b^n, 

r Ai (2.12) 

9aa f,, 

where the index a runs over the set of particles localized in W space, except for n. 



3 A CSW-like Expansion for Gravity 

In [3], a derivation of the CSW expansion for the A/" = 4 super Yang-Mills NMHV amplitudes from 
the Grassmannian is presented using the "relaxing (5-fuctions" procedure described in section 3 of [3] . 
Here we apply the same procedure to gravity, but instead of working in momentum-twistor space, we 
start from the link representation formula (2.8) written in momentum space. 

Let us here summarize the procedure. To start with, it is useful to rewrite the products in (2.8) 

as 



^12 '~-la'~-2aWn—lWn — \^12 ) 



d-1 / \ d-1 

CiaC2a I ( Cm—lCrn \ i 



SO that 



7-)^^— 1 n 
^12 



{cin-lClnC2n-lC2nY 



^12 ClaC2aCrn—lCr 



rj^l,2 
aj^n—l,n 



d-1 



d-1 



J\ ClaC2a Cm-lCr 



yaf^n.n—l 



Then, we can express (2.8) as 



Mn,d = / n 7^ n ('^i^'^aa 

r,a a^n,n-l 



d 2 d 2 (-^12""'^") 



r=/=l,2 



{Cln- 1 Cl„C2n- 1 C2n )^ 



a^n,n— 1 



(3.1) 



For the NMHV sector, we have d = 2 and d = n — A; so the number of integration variables is 
(d — l)(d — 1) = n — 5. Also, we choose the gauge 



'10 Ci4 Ci5 ... Cln 
1 C24 C25 ••• C2„ 
.001 C34 C35 ... C3„ , 



(3.2) 
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such that r, s = 1, 2, 3 and a, 6 = 4, 5, . . . , n. Then, (3.1) written in momentum space reduces to 

M... . /n f^,.„-.30"-^J^££;^.- (^) (M) 

(3.3) 

n-2 1 3 n ^ ^ 

X n T^i2^^ n -^'(^^ + ^™^") n - -^'^™) • 

a=4 ''^ri-lna a=4 

Now, starting from (3.3), we perform the fohowing steps: 



1. Pull out momentum conservation 

3 



n 5\K + CraX") = (23)^5^ A,,V ; (3.4) 

r=2 \i=l / 

2. Split the remaining Ai "two component" (5-functions into two "one-component" ^-functions by 
projecting it on two arbitrary linearly independent spinors [X\ and 

5''{\i + ciaX") = [XY]5{[IX] + cia[aX])5{[lY] + cia[aY])- 

3. Solve the system with all A (^-functions together with (5([iy] + cia[ay]), and relax 5([1X] + 
cia[aX]), that is to replace the 5-function by one over its argument and treat the integral as a 
contour integral around the point where this argument is zero: 

5{[lX]+cia[aX])- ^ 



[lX] + cia[aXy 

4. Since we relaxed one 5- function, we increase by one the number of integration variables ti, i = 
1, 2, . . . , n — 4. The amplitude M.n,2 is obtained by carrying the integrations on a contour defined 
around the zeros of the n — 5 Veronese polynomials plus the zero of [IX] + cia[aX]. However, 
relaxing the (5-function means treating this point as a pole, so we can use the global residue 
theorem to "blow up the residue" and write the amplitude as minus the sum of residues of the 
other poles of the integrand. Being more precise, suppose we perform n — 5 integrations to solve 
all Veronese constraints. Then, we are left with one integration, say in the variable ti, and the 
integrand is a function of ti. To obtain the amplitude, we should compute the residue due to the 
only remaining constraint [IX] + cia(ti)[aX] = 0. However, if the integrand contains m other 
poles for t^ j, j = 1, . . . ,m, then 



M.n,2 = — Res [Integrand] 



(3.5) 



where these poles correspond to 3 x 3 minors in the denominator. 

These steps are the core of the CSW expansion for super Yang-Mills amplitudes, and in this work we 
apply the same procedure for the gravity formula (3.3) to calculate an analogous expansion for 6- and 
7-graviton NMHV amplitudes. 

Let us now explain how this expansion (3.5) coincides with CSW. We discuss the localization in 
the space of 3-vectors, C-space, which are the columns of a 3 x n matrix built from three n-vectors 
(3.2). Both Yang-Mills and gravity amplitudes in the NMHV sector localize on pairs of lines in the 
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C-space. For Yang-Mills, this localization is equivalent to the localization in twistor space, while for 
gravity the situation is different, as we know even the MHV amplitudes in gravity are not represented 
by lines in twistor space. 

Veronese map is a one from G(2, n) to G{k, n) which transforms world-sheet variables to fe-vectors, 
columns of a A; x n matrix. Veronese polynomial, defined in (2.10), when is set to zero, is a constraint 
for the columns of the C-matrix to be the same as columns of this k x n matrix. 

In the computations of residues of (3.3), the C-matrix is evaluated on each of the poles of the 
integrand. This 3 x n matrix, after a rescaling of one of its rows, can be depicted on a 2-plane. Since 
we have n — 5 Veronese polynomials in (3.3), not all but some of the solutions make the points of 
the rescaled C-matrix lie on a conic. There also exist some non-conic configurations. The numerator, 
however, vanishes on the solutions leading to these configurations, so projects them out. Moreover, 
at any time, as there is a minor from the denominator set to zero, three of the points have to be 
collinear. In general, any five points define a unique conic. When three of these points are collinear, 
the conic is reducible and the configuration will be a degenerate conic which has three points on 
one line and the others spreading on any of the two lines. So, in summary, on the support of n — 5 
Veronese polynomials and the numerator we only see conic configurations. Each of the minors in the 
denominator then tells us which 3 points must be collinear. The outcomes are the terms in (3.5) 
which we define to be CSW terms. 

Here we refer to our expansion not as a "CSW expansion for gravity" , but as a CSW-like expansion. 
The reason is that "the canonical" way of generating the MHV expansion for gravity amplitudes is 
through Risager's procedure. It is not obvious a priori, however, that the procedure shown here will 
correspond to Risager's as happens for super Yang-Mills. That point being clear, we will refer to the 
channels as CSW-terms rather than CSW-like-terms for convenience of the notation. 



4 6-Graviton Computation 

For the 6-graviton NMHV amplitude, we substitute n = 6, d = d = 2m (2.8). As already said, 
we choose a gauge fixing such that the indices r, a take values r = 1,2,3 and a = 4,5,6. For 
computational purposes, here we write the amplitude not in twistor space but in momentum space, 
so the factor exp {icra^^y^"') from (2.8) is Fourier transformed into a product of J-functions. Then 
(2.8) is written as 



^^(l-,2-.3-.4^5^6^) = JuJ-^(H!l)V^' (^) (g) 



r=4,5,6 



6 3 



(4.1) 



V 



a=4 



r=l 



where 



^ c^^c^e^c^ and V^VfH 



(4.2) 

C15C26C16C25 

The functions <f>^'^^ entering the formula (4.1) are determinants of arbitrary 2x2 minors of the following 
rank-2 matrices: 



(12> 


(13> 


(12> 




{13} 




m 


(12> 


(23> 


Ff 

(23) 


rr56 




rr56 

12, > 


mi 


, mi 


(13) 




(23> 




(13> 


HE 




mi 




rr56 
^13 



(4.3) 
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There are 9 integration variables and 12 bosonic 5- functions, four of which provide momentum- 
conservation (3.4). After pulhng out the momentum conservation the integral (4.1) contains 8 con- 
straints and 9 integration variables, so it should be complemented by the choice of an integration 
contour. The amplitude is given by integration along the contour encircling V = 0. 

Now let us start the transformations leading from the integral representation (4.1) to the CSW-like 
expansion of the gravity amplitude. First, according to the procedure explained in section 3, we split 



^^(Ai + C14A4 + C15A5 + cieAe) = [XY]5 {[IX] + ci4[4X] + ci5[5X] + cie[QX]) 

x6{[lY]+cu [4r] + ci5 [5Y] + cig [6^] ) , 

and relax 5 ([IX] -|-cia[aX]). All the factors in the denominator can then be combined into two 
functions: 

fi = ci5Ci6C25C26C34(156)(256)(345)(364) {[IX] + cia[aX]) , (4.6) 
/2 = V. (4.7) 

The integral is now in two complex variables and the amplitude is then given by a multi-dimensional 
residue at the poles where both [IX] + cia[aX] and V are zero. 

According to the global residue theorem, we can also write the amplitude as minus the sum of 
all other residues in which V and one of the factors in /i other than [IX] -|- cia[aX] are set to zero. 
Explicitly speaking, denoting a residue by the factors set to zero to compute it, then 

Me,2 = {V, [IX] + cia[aX]} = -{V, C15} - {V, cie} - {V, C25} - {V, cse} - {V, C34} 

-{V, (156)} - {V, (256)} - {V, (345)} - {V, (364)}. 

Notice that, since V is a polynomial of degree 4, each term splits yet into 4 terms. Notice also that 
many residues contribute to the same configuration, e.g., {(156), (126)} and {(256), (126)} shown in 
figure 2. The way to carefully deal with this and take care of all contributions for a given diagram 
will be explained in section 4.2. 




Figure 2: A configuration associated with {(156), (126)} and {(256), (126)}. 



Carrying on the computation, the (5-functions which we keep, impose the following set of seven 
equations 

' [ly] + ci4[4y] + ci5[5y] + ci6[6y] = o, 
^ A4 = C14A1 -I- C24A2 + C34A3, 
A5 = C15A1 + C25A2 + C35A3, 

^ Ae = cigAi -I- C26A2 + C36A3, 
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or equivalently 

' [lY] + cu[4Y] + ci5[5Y] + ci6[6y] = 0, 
ci4(12)+C34(32)-(42) = 0, 

C24(21)+C34(31)-(41) =0, 

< ci5(15) + C25(25) + C35(35) = 0, (4.10) 

C25(21)+C35(31)-(51) =0, 
C16(16)+C26(26)+C36(36) = 0, 
.C26(21)+C36(31)-(61)=0. 

The last transformation is performed for computational convenience and multiplies the Jacobian by 
a factor of (21) (51) (61). 

Let us denote by S the system of equations formed by (4.10) and two additional equations: V=0 
and (abc) = 0, where (abc) is one of the minors in the denominator. Then, an individual residue is 
given by the following expression evaluated on the solution of the system S: 

1 ([45][56]ci5C25(364) + cyclic(4, 5, 6)) ((12)(23)c25C26(256) + cyclic(l, 2, 3)) jabc) [XY] 
J (135)(136)(235)(236)(124)(156)(256)(345)(364) {[IX] + cia[aX])' 

(4.11) 

where J is the Jacobian of solving the system S. The precise way of calculating the residues will be 
explained later in section 4.2. 



4.1 Analytical Computation of {(135), (246)} 

We calculate the residue of the integrand at the point {(135), (246)}. To see that this is indeed a simple 
multi-dimensional pole we represent the Veronese polynomials in the form V = (123)(345)(561)(246) — 
(234)(456)(612)(351). Setting (135) to zero makes /i vanish. Setting additionally (246) to zero sets 
both /2 and V to zero and, as one can check, does not add any zeros among the factors in the 
denominator. 

Using Mathematica to find the solution of the system (4.9), substituting it into (4.11) and simpli- 
fying the result, we obtain 

[46](13)6(2|4 + 6|5](2|4 + 6|y]5 

(15)(24)(26)(35)(46)(l|3 + 5|y](4|2 + 6|y](6|2 + 4|y](5|l + 3|y](3|l + 5|y](pi+p3+p5)2' ' 

One can see that this result coincides with a Risager's term which has the form 

A^L(i-,3-,5+,P+)-^A^jj(2-,4+,6+,-P-), (4.13) 

where P = P2 + Pa+ Pe is the momentum flowing through the link which is set on-shell and P is the 
Risager-deformed momentum P = P2 + P4+ Pe as in (1.2). 



4.2 Numerical Check of All Residues for 6 Gravitons: Taming Sin- 
gular Configurations in Multi-Dimensional Residues 

In order to numerically calculate the other channels, an immediate application of the global residue 
theorem leads to a computational problem: the part of the integrand besides the two zero factors in 
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the denominator, which is supposed to be finite, can become ambiguous. This happens because the 
numerator and the denominator go to zero simultaneously and do not cancel with each other. 

Consider for example the 2-4 channel shown in figure 3 in which particles 1 and 4 belong to Ml 
and particles 2, 3, 5 and 6 belong to M.r- 




Figure 3: An example of a 2-4 channel. 

One way to generate this configuration is by setting V and (256) to zero. Looking at the denominator 
of (4.11), however, we see that (235) and (236) are "accidentally" set to zero as well. For the residue 
to be finite, the numerator of (4.11) should also go to zero at the solution. Restricting ourselves 
to the submanifold of G{3, 6) defined by V = equation, we can write all minors that go to zero 
as a multiple of say (235). More precisely, in this example we can use V = (462) (235) (514) (631) — 
(623) (351) (146) (254) = 0^ to write (236) oc (235). Playing this game with all the minors that go to 
zero in the integrand, we note that both the numerator and denominator vanish as (235)^, so they 
cancel and we can finally obtain a finite residue. 

Although this method gives the correct results in this example, it requires an individual treatment 
of each residue, making the computations rather laborious. More importantly, it cannot necessarily 
be applied to higher-point computations. For this reason, we use a procedure which can be directly 
applied in all cases which consists of the use of a regulator. 

The singular residues are points in the Grassmannian on which the denominator of the integrand 
has a zero of second order or higher. The idea is to add to each minor in the denominator a small 
constant e, the regulator, 

{d^ ^ (a6c) + e' ^^'^^^ 

which separates the higher order zeros in many zeros of first order. Also to make sure that the regulated 
minors do not accidentally vanish simultaneously, we may add different values of e to different minors. 
After this, all poles are simple and the global residue theorem can be applied without restriction. 
Moreover, there is no need to regulate the Veronese polynomials, since they do not factorize when 
minors are regulated (e 7^ 0). 

The Veronese polynomials and each of the regulated minors in the denominator make a system 
of equations set to zero. On each solution of each system of equations, there might be some minors 
which nearly vanish. We compute all possible minors on each solution and collect the almost-zero 
ones. These minors now define the localization of particles on the two crossing lines which makes a 
CSW channel. As an example, in figure 3, the only minors which are close to zero are (235), (236), 
(256) and (356). Hence, particles 2, 3, 5 and 6 lie on a line. And, there is also a line passing through 
any two points; 1 and 4 here. As it is also clear from this example, many solutions may contribute 
to a given localization. Here, residues of the integrand on the solutions of the following systems of 
equations contribute to this channel: V = (235) + e = 0, V = (236) + e = and V = (256) + e = 0. In 
the end, all contributions must be summed to give the corresponding CSW term. 

"'^Keep in mind that the condition V = is permutation invariant, so we can make a choice of ordering of particles which 
is the most convenient for our purposes. 
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This way, we computed the numerical values of the 21 configurations of 6 gravitons, depicted in 
figure 4. The numerical values obtained are compared to the known corresponding Risager terms. As 




Figure 4: 6-graviton channels 



a result of this comparison we found, as happens to Yang-Mills, a complete agreement between our 
procedure and Risager's. Numerical results are given in the appendix for both 6- and 7-graviton cases. 
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5 7-Graviton Computation 

For n = 7, d = 2 and d = 3, so the general formula (2.8) is 
A^(l-,2-,3-,4+,5+,6+,7+)= / J] 



dCra (367)^ 

C36C37 1 



_, „ „ Cra (C16C17C26C27 
a=4,.,7 



(5.1) 



\J^12/ \^rs / 1^467 "^567 „ 

Using the definitions of (j), one can rewrite the integrand as a rational expression whose denominator 
is 



(136) (137) (236) (237) (124) (125) (167) (267) (345) (346) (347) (356) (357) V]^^ Vg^gf • (5-2) 

According to the general procedure of section 3, we relax one component of the Ai 5-functions. After 
this, we get an additional factor of [IX] + X^aCia[aX] in the denominator. To apply a global residue 
theorem these factors should be split into three functions: 

/i = (136)(137)(236)(237)(124)(125)(167)(267)(345)(346)(347)(356)(357) ( [IX] + ^ cia[aX] ) , 



^123 



/3 — 



467) 

il23 
567- 



(5.3) 



The amplitude is the residue of the integrand at the point where /2, /s and the last factor in /i go to 
zero. Like in the 6-particle case, we get different CSW terms by computing residues at points where 
/2, /s and one of the minors in /i go to zero. The GRT then tells us that all the CSW terms sum to 
the full amplitude, as it naturally should be: 



M 



7,2 



[IX] + cia[aX]j , /2, /aj = -{(136), /2, h} - {(137), /2, /a} - • • • - {(357), /a, fs}. 



(5.4) 



Notice that every term in the right-hand side can in principle contain more than one CSW term, 
because the system 



minor = 0, 

Vill = 0, (5.5) 

;123 
567 



i;123 _ r. 



can have more than one solution. 

To compute the residues we perform a procedure similar to the one described in section 4.2. First 
we introduce the regulator to every minor in (5.2) and write a system of linear equations composed 
of the (5-functions and three equations setting a minor and the two V's to zero. After this, all residues 
are free from any accidental ambiguity and out of the 13 x 4 x 4 solutions, many of them sum to give 
the 45 possible configurations. 

If the configuration is such that these two Veronese polynomials are equal to zero, but some other 
Veronese polynomials which can be made of particles 1, ... ,7 are not zero, then this configuration 
does not contribute, since the corresponding residue is zero (such configurations are called spurious 
solutions). 
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6 The General n NMHV Case: Power Counting 



All the reasoning presented in section 3 relies on the assumption that the global residue theorem is 
applicable and our aim now is to verify when this is the case for gravity. To do so, we analyse how 
the integrand behaves for large values of the integration variables, that is, we look for n such that 
there is no pole at infinity. 

In order to perform the power counting, consider formula (3.3) which we repeat here for conve- 
nience: 

U - / TT ^^^'^ r\(r n (-^"2'^")" ^ ,(2) ( A ,(n-4) (\(^^X\ 



r,a — a=4 y-m- i'-ln'^zn- i'^zn J V-^rs / \-'-'12. 

n-2 

X 



n— 2 ^3 n 

n n + ^-^") n - 

a=4 '^n-lna ^=1 a=4 



We know from section 2.4 of [18], that each minor is at most of degree one. We also know that in the 
link representation we can write each link variable and the factors D^2^^ minors, so they are also 
linear in ti, and keep in mind that each Veronese polynomial is of degree 4 in tj. 



Then, let us consider one variable t and look at how each factor in the integrand of (3.3) scales 
with it: 



j.9-3n 
~ t , 



\n-6 .2n-12 



^ \n—b J.Z 
(CSn-lCSn) ~E 

+71 — 4 

^n-12 



(ci„-lCi„C2n-lC2n)^ 
, (n-4) ( M\ _ ,3(n-4) 

\HfJ 

n i ~ t-'"-'>. 

Relaxed (5- function ~ 

Summing up, we get that the integrand scales as The condition for the global residue theorem 

to apply is that the integrand should decay faster than for each integration variable, that is 

—n — 2 ^ —n + 4. This condition is satisfied Vn, so we conclude that the CSW-like expansion will be 
valid for all NMHV graviton amplitudes. 



7 Conclusions 

Building up on the recent results for A/" = 8 supergravity amplitudes — Hodge's MHV formula, the 
Cachazo-Skinner formula and its link representation — we carried on the appliance of the well-known 
techniques used for super Yang-Mills amplitudes to explore the new gravity territory. 

From the link representation, we derived a CSW-like expansion for pure 6- and 7-graviton NMHV 
amplitudes. It was not obvious, however, that this expansion would coincide with Risager's, the 
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known way up to now to produce an MHV-vertex expansion. We observed a complete agreement 
between the two expansions in the cases we worked on. 

The most exciting feature of the new expansion is that in principle it works for any n for d = 2, 
while Risager's method fails for n > 12. Further studies are necessary to verify if and when the two 
expansions stop agreeing. One possible reason for the disagreement of the two expansions is CSW-like 
configurations which are not the product of two MHV amplitudes but are still non-zero residues of 
the link formula. In other words, we may have a pair of lines in the C-space with all negative helicity 
particles on one line. This is obviously not a Risager's term. For small n, as was said, the numerator 
does not allow such configurations in the expansion; in fact, the residue will be zero. But, for higher 
n one can still have a valid expansion of residues, as the power counting shows, but there exist more 
terms which are not MHV x MHV. These extra terms are in fact residues at infinity of the Risager's 
method. It would be interesting to explore these extra terms more precisely. 
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A Numerical Values for the 6- and 7-Graviton Residues 
and Risager Terms 

Here we present the numerical values for different channels given by the Risager diagrams and our 
CSW-like terms for a given set of external data for the 6- and 7-graviton cases. As described in section 
4.2, the residue computations are based on the use of a regulator e for the minors of the integrand. 
We chose e to take the value 0.0000001. Also, note that the amplitude is the sum of Risager terms but 
is minus the sum of residues as is known from the residue theorem. So, here, signs of the values of the 
two columns are opposite. The biggest errors in the 6- and 7-graviton computations are 1.6 x 10~^% 
and 1.4 X 10~^% respectively. 



Our data for the 6-graviton computation are as follows: 
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Table 1: Numerical Values for 6-Graviton Computation 



16 



(^Vl o-n-npl 


± LlOCl&,d -LC/l 111 


R pQi H n p 


Error 


|i,ZJ-|J,4,0,D, 1 J- 


— 1 4'?n9 V 1 (\~'^ 

— 1.4:OU^ A lU 


1 4^97 V 1 O"'^ 


1 7 
1. ( 


X 


lU 


A) 




—0 nni 9fi^nfi4fi 


001 9f\^f)f\9^ 


9 n 
z.u 


X 


1 

lU 


/o 


/I 4IJ9 fi 71 


4 075 X 1 n^^ 


—4 099 X 1 0^^ 


P. 1 
0.1 


X 


lU 


/o 






—0 0000^49(^44^"^ 


1 8 
1.(5 


X 


1 C\~ 
lU 


/o 




\j .yjyjyjO'^o j-OijO 


000^^451 ^Q^^ 


0.0 


X 


lU 


/o 


/I 7119 4 ^ fil 


—0 DDI 1 R^^^zL9 


001 1 J^"^^^^)^ 
u.uu 1 1000000 


1 /I 
1.4 


X 


lU 


/o 


■[9 '^in 4 ^ (S 71 


— n 00941 Ql 9(S1 


00941 Ql 9'^9 


1 9 
l.Z 


X 


lU 


/o 


19 4li1 ^ F, 7\ 


9 R(^Q^70 V 1 0""^ 

Zi.OUc/tJiU A ±U 


—9 R(SQ^^^ V 1 O"'^ 

Zi.OUc/OOO A lU 


9 
o.z 


X 


1 

lU 


/O 


19 ^in 3 4 fi 71 




— 9 ^RQ X 10"'^ 


1 /I 
1.4 


X 


1 n- 

lU 


/O 


{9 fil-Tl 4 ^ 71 


— n nnni "^si n87 


n nnni "^si ns^ 


1 /I 

1.4 


X 


1 n- 

lU 


/o 


/9 71/1 4 fil 


n 091 81 9 
u. uz 10 IZ 


— n n9i snn 


1.0 


X 


lU 


/o 


i'? 4li1 9 (S 71 


—0 nnnn'?sqsfi4 

u.uuuuooyou4 


n nnnn'?SQsn7 


7 7 
( . ( 


X 


lU 


/o 


i'? ^lil 9 4 (S 71 


n nnnn4f^447 


— n nnnn4^4^n 

U.UUUU4tJ4tJU 


0.0 


X 


1 n- 

lU 


/o 


■T'? (Sin 9 4 71 
0111,^,4,0, r ]■ 


— n nn9i n89'?7 


n nn9i n894'? 


9 8 
Z.o 


X 


1 n- 

lU 


/O 


-T'^ 71/1 9 4 1^ fil 


—0 nn'?9fi8'i8 


n nn'?9fi8'i4 

U.UU0^U004 


1 9 
l.Z 


X 


1 n- 

lU 


/o 


/I 9 41/'? ^ fi 71 


— Q "ii 1 Qfi V 1 


Q "ii 141 V 1 n~^ 


0.0 


X 


1 n- 

lU 


/o 


-Tl 9 '^11'? 4 71 


— '? '?fi'^'?i 1 V 1 

O.OUOOll y\ lU 


'? '?fi'^'?fi9 X 1 n~8 


1.0 


X 


lU 


/o 


/I 9 fill'? 4 ^ 71 


o.y/zyou a lu 


— ^ Q79Q44 X 1 n^8 

— 0.y(Zy44 A lU 


l.U 


X 


lU 


/o 


-Tl 9 711'? 4 ^ (Si 


— n9n48Sfi 


n n9n48^7 


1 /I 
1.4 


X 


1 n- 

lU 


-40;/ 
/o 


/I ■? 41/9 fi 71 


— u.uuuoouuyyi^iioo 


n nnn'?'?nnQQi4i '?9 

U.UUU00UUC/C/1410Z 


o.U 


X 


lU 


/o 


/I ■? "11/9 4 fi 71 




n nnin8^Q949 


1 n 

l.U 


X 


1 n- 
lu 


/o 


/I '? (Sl/9 4 ^ 71 


— n nn^8nn'?n9n 


n nn^snn'^nfi'^ 


7 8 
( .0 


X 


1 n- 

lU 


-70/ 

/o 


/I '? 71/9 4 ^ (Si 


n nnn7fi7i4i 

U.UUU 1 1 141 


— n nnn7fi7i '?q 

— U.UUU ( ( loy 


P. 

z.o 


X 


lU 


-60/ 

/o 


/I 4 ^1/9 '? (S 71 


1^ 89'?'?Q'?7Q V in~6 

O.OZOOc/O/c/ A lU 


— ^ 89'?'?Q'?f;8 X in^6 

— O.OZOOc/OOO A lU 


0.0 


X 


lU 


/o 


/I 4 (Sl/9 '? 71 


— n nnnn4Sfinsi 1 

— U.UUUU40UU01 10 


n nnnn4Sfinsi 91 


1 
l.Z 


X 


lU 


/o 


/I 4 71/9 "? fil 


— n nni fin^i 9nn 


n nni fin"^! 997 


1 7 
1. / 


X 


1 

lU 


/o 


/I 1^ (Sl/9 '? 4 71 


n nnni 4RfiRn4 


— n nnni 48fi8ni 

U.UUUI 40 UOU 1 


n 
z.u 


X 


lU 


-60/ 

/o 


/I f> 71/9 '? 4 (Si 


n nn9'?fi^ss4fi 

U.UUZ0000040 


— n nn9'?fipi88fin 

— u.uuzouooouu 


o.y 


X 


lU 


-70/ 

/o 


/I (S 71/9 '? 4 "^1 


n nnni snnfini 

U.UUUlOc/UOUl 


— n nnni snnfinn 

— U.UUUlOc/UOUU 


0.0 


X 


1 0" 
lU 


/O 


/9 '? 41/1 ^ fi 71 

■[Z,0,4J-^[l,0,0, ( } 


n nn45n59957 


— n nn45n59935 


4 Q 
4.y 


A 


1 0" 

lU 


/o 


So Q c;l/i 4 fi 71 


— n nnnnQ9fififiQ 


n nnnnQ9fifiQfi 


9 Q 

z.y 


X 


1 n- 

lU 


-50/ 

/o 


/9 '? fil/1 4 ^ 71 
■[z,o,0|^[i.4.o, / 1 


0099'^'^?'^^ 


—0 0099'^'^7fi0 

U.UUZ, ZiOO 1 uu 


z.z 


X 


lU 


-60/ 

/o 


/9 '? 71/1 4 ^ fil 




—4 941 ^491 V 1 0~6 


8 

0.0 


X 


lU 


-80/ 

/o 


/9 4 f>l/1 '? fi 71 


—A 00R10Q7 V 10^9 


4 00^1 0^1 V 1 0^9 


4 n 

4.U 


A 


1 n^ 

lU 


/o 


/2 4 fil/1 '? 5 71 

■[Z,,4,UJ-^[l,0,0, ( J- 


—0 0007^^8(1578 


n nnn738fi558 


9 7 
Z. / 


X 


1 n- 

lU 


/O 


/9 4 71/1 "? fil 


008^41400 


— n nn8^4i499 


1.0 


X 


lU 


-60/ 

/o 


/9 ^ fil/1 '? 4 71 
■[Z,0,0J-^[l,o,4, / 1 


0001 ^^^^1 ^ 
u.uuui 0000 10 


—0 0001 ^^^^1 ^ 
— U.UUU 10000 10 


1 fi 
1.0 


X 


1 n- 
lU 


-60/ 

/o 


/9 ^ 71/1 '? 4 fil 


—0 00"^! Q^(^Q 

— U.UUOlolc/OUc/ 


OO'^l Rl Q^(S1 

U.UUOlOlc/OUl 


2.5 


X 


10" 




■[2,6,7111,3,4,51 


-0.0005354038 


0.0005354037 


1.9 


X 


10" 




{3,4,5}{1,2,6,7} 


0.000044443578 


-0.000044443556 


5.0 


X 


10" 




{3,4,6}{1,2,5,7} 


0.002633 


-0.002671 


1.4 


X 


10- 




{3,4,7}{1,2,5,6} 


0.00395962 


-0.00395958 


1.0 


X 


10- 




{3,5,6}{1,2,4,7} 


-0.006673719 


0.00667372 


1.2 


X 


10- 




{3,5,7}{1,2,4,6} 


2.363790 X 10"^ 


-2.363757 x 10"^ 


1.4 


X 


10- 




{3,6,7}{1,2,4,6} 


-0.000026712880 


0.000026712871 


3.4 


X 


10- 





Table 2: Numerical Values 



ir 7-Graviton Computation 



References 

[1] S. He, "A Link Representation for Gravity Amplitudes," arXiv: 1207.4064 [hep-th] . 

[2] F. Cachazo, L. Mason, and D. Skinner, "Gravity in Twistor Space and its Grassmannian 
Formulation," arXiv: 1207.4712 [hep-th]. 

[3] N. Arkani-Hamed, J. Bourjaily, F. Cachazo, and J. Trnka, "Local Spacetime Physics from the 
Grassmannian," JHEP 1101 (2011) 108, arXiv:0912.3249 [hep-th]. 

[4] K. Risager, "A Direct proof of the CSW rules," JHEP 0512 (2005) 003, 
arXiv : hep-th/0508206 [hep-th] . 

[5] F. Cachazo, P. Svrcek, and E. Witten, "MHV vertices and tree amplitudes in gauge theory," 
JHEP 0409 (2004) 006, arXiv:hep-th/0403047 [hep-th]. 

[6] R. Britto, F. Cachazo, and B. Feng, "New recursion relations for tree amplitudes of gluons," 
Nucl.Phys. B715 (2005) 499-522, arXiv:hep-th/0412308 [hep-th]. 

[7] R. Britto, F. Cachazo, B. Feng, and E. Witten, "Direct proof of tree-level recursion relation in 
Yang-Mills theory," Phys. Rev. Lett. 94 (2005) 181602, arXiv:hep-th/0501052 [hep-th]. 

[8] M. Bianchi, H. Elvang, and D. Z. Freedman, "Generating Tree Amplitudes in N=4 SYM and N 
= 8 SG," JHEP 0809 (2008) 063, arXiv : 0805 . 0757 [hep-th]. 

[9] P. Benincasa, C. Boucher- Veronneau, and F. Cachazo, "Taming Tree Amplitudes In General 
Relativity," JHEP 0711 (2007) 057, arXiv:hep-th/0702032v2 [HEP-TH]. 

[10] E. Conde and S. Rajabi, "The Anomaly of the Twelve-Graviton Next-to-MHV Risager 
Amplitude," JHEP 1209 (2012) 120, arXiv: 1205 . 3500 [hep-th]. 

[11] E. Witten, "Perturbative gauge theory as a string theory in twistor space," 

Commun.Math.Phys. 252 (2004) 189-258, arXiv:hep-th/0312171 [hep-th]. 

[12] A. Hodges, "A simple formula for gravitational MHV amplitudes," arXiv: 1204. 1930 
[hep-th] . 

[13] F. Cachazo and D. Skinner, "Gravity from Rational Curves," arXiv: 1207.0741 [hep-th]. 

[14] M. Bullimore, "New Formulae for Gravity Amplitudes: Parity Invariance and Soft Limits," 
arXiv: 1207.3940 [hep-th]. 

[15] N. Arkani-Hamed and J. Kaplan, "On Tree Amplitudes in Gauge Theory and Gravity," JHEP 
0804 (2008) 076, arXiv : 0801 . 2385 [hep-th]. 

[16] P. Griffiths and J. Harris, Principles of Algebraic Geometry. Wiley Classics Library Wiley 
Classics Library. Wiley, 2011. 

[17] R. Roiban, M. Spradlin, and A. Volovich, "On the tree level S matrix of Yang-Mills theory," 
Phys. Rev. D70 (2004) 026009, arXiv:hep-th/0403190 [hep-th]. 

[18] N. Arkani-Hamed, F. Cachazo, C. Cheung, and J. Kaplan, "A Duality For The S Matrix," 
JHEP 1003 (2010) 020, arXiv: 0907. 5418 [hep-th]. 77 pages, 19 figures. 



18 



